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Abstract. Detailed first estimates are presented of an- 
gular broadening in the outer heliosphere due to scat- 
tering of radio waves by density irregularities. The ap- 
plication is to the 2-3 kHz radiation observed by Voy- 
ager. Two plausible turbulence models, which account 
very well for scattering within 1 AU, are extrapolated 
beyond 10 AU. Both models predict significant angular 
broadening in the outer heliosphere, accounting semi- 
quantitatively alone for the source sizes inferred from 
roll modulation data. Predictions are presented for ra- 
dial variations in the apparent source size if scattering 
is important. Comparisons with available data argue 
that scattering is important (and indeed is the domi- 
nant contributor to the apparent source size) and that 
the radiation source is located in the outer heliosphere. 
Other evidence that scattering is important, such as the 
fluctuations in apparent source direction and intensity, 
are also identified. The effects of scattering should be 
included in future analyses of the 2-3 kHz emissions. 

Introduction 

Density irregularities are well known to scatter radio 
waves in solar system and astrophysical plasmas, affect- 
ing the apparent angular size, intensity, and direction 
of the radio source, as well as the temporal and spectral 
variability of the radiation [e.g., Hollweg , 1970; Stein- 
berg et ai, 1971, 1985, 1989; Rickett, 1977]. Radiation 
generated or propagating close to the electron plasma 
frequency f p is strongly affected by scattering. Indeed, 
analyses of all known emissions generated near f p and 
2/p in the corona and the solar wind within 1 AU show 
that scattering is very important in interpreting the an- 
gular size, directivity and brightness temperature of the 
radiation. These sources are type III solar radio bursts 
in the corona and solar wind [Steinberg et aL, 1971; 
1985] and 2 f p radiation generated in Earth’s foreshock 
[Lacombe et aL, 1988]. Scattering also significantly af- 
fects AKR observed in the solar wind [Steinberg et aL, 
1989] even though it is generated and propagates well 
above f p from its auroral sources. 

The Voyager spacecraft have observed radiation prop- 
agating near 2-3.5 kHz in the outer heliosphere [Kurth 
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et aL, 1984; Kurth, 1988; Gurnett et aL, 1993; Gumett 
and Kurth, 1994] at frequencies ~ 1 - 10 times the local 
solar wind plasma frequency and ~ 1 - 5 times likely 
plasma frequencies in the heliosheath. Existing theories 
all predict generation of the radiation near f p or 2/p in 
differing outer heliospheric sources [e.g., Cairns et aL, 
1992; Gurnett et aL, 1993]. Modulations of the radia- 
tion intensity near 3 kHz during Voyager roll manoeu- 
vers are observed and interpreted in terms of the appar- 
ent size of and direction to the radiation source [Kurth, 
1988; Gurnett et aL, 1993]: near 48 AU these data are 
consistent with a maximum half-angle ~ 60° ~ 1 rad 
for the source. It is relevant to note the absence of 
modulations near 2 kHz and the significant variation 
in apparent source direction [Gumett et aL, 1993]. In 
addition, Cairns et aL [1992] ascribed fluctuations in 
the radiation’s intensity and frequency fine structure to 
scattering and/or ducting. Based upon these properties 
and the analogies with inner heliospheric radio sources, 
it is very plausible that scattering off density irregular- 
ities might be important in determining the properties 
of the Voyager 2-3 kHz radiation. This paper provides a 
first, detailed, semi-quantitative investigation into how 
scattering affects the source’s apparent angular size. 

Scattering rates near 50 AU: first 
calculations 

Two formalisms exist to calculate the angular broad- 
ening of sources due to scattering: a geometric optics 
approach that has been used in most solar system ap- 
plications [e.g., Hollweg, 1970; Steinberg et aL, 1971, 
1985, 1989; Lacombe et aL, 1988], and the approach 
based on the parabolic wave equation (hereafter PWE) 
that is favored in astrophysical applications [e.g., Lee 
and Jokipii , 1975; Rickett, 1977; Spangler and Arm- 
strong, 1990]. The physical difference is that diffraction 
is not included in geometric optics. In both cases the 
apparent half- angle of the source, $appi I s 

0 a pp « yj 0,n + ^sc (!) 

Here, 0 m and 0 SC are half-angles for the intrinsic source 
size and the broadening contribution due to refrac- 
tive scattering by density irregularities, respectively. 
In both approaches 0 K depends on a path integral 
dz from the source to the observer of an integrand 
that depends on the radiation frequency /, the local 
(steady-state) plasma frequency f p 0 (z), and the three- 
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dimensional wavevector spectrum of density turbulence 
P n (k, z). Further shared assumptions are that the in- 
ner scale /[ = q ~ 1 of the turbulence greatly exceeds the 
radiation’s free-space wavelength A (= c/f) and that 
the outer scale l Q is much less than the path length L. 

Geometric optics is conventionally a viable approxi- 
mation when [e.g., Tatarskiy 1971] 

L Z if/ A . (2) 

Taking the thermal proton gyroradius as the inner scale 
/in and radiation at 2-3 kHz, near 50 AU /i n ~ 10 6 m, 
and so by (2) geometric optics is viable only for path 
lengths L £ 10 7 m ~ 10“ 4 AU. However, the minimum 
plausible distance to the source is ~ 50 AU since the 
Voyager spacecraft are separated by a comparable dis- 
tance but observe very similar radiation. Geometric op- 
tics therefore appears not viable in the outer heliosphere 
[J. W. Armstrong , personal communication, 1995]. Con- 
sequently the PWE formalism is used in this paper. 

The turbulence spectrum P„ (k) is currently unknown 
experimentally or theoretically beyond ~ 5 AU. Unan- 
swered questions include how well WKB theory applies 
and how MHD turbulence driven by interstellar pick- 
up ions affects the spectral index and level of P n (k) 
at small scales (where scattering is most important). 
It is therefore necessary to introduce plausible semi- 
quantitative models in order to make first estimates 
of whether scattering is important in the outer helio- 
sphere. These models must account adequately for ob- 
servations within 1 AU and be extrapolated plausibly. 

For a Gaussian distribution of density irregularities, 
defined by a characteristic scale h and a fluctuation 
level An/n, and generalizing the PWE results of Lee 
and Jokipii [1975] and Rickett [1977] to include finite 
fpo(z)/ fy it may be shown that 


e 


2 

sc 
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Scattering thus increases for smaller h, and larger ratios 
fp/f and An/n. Taking h equal to the thermal proton 
gyroradius r p ~ 100 km and An/n = 10~ 2 , choices 
that are both conventional and plausible at 1 AU [e.g., 
Hollweg, 1970], and writing 0 2 c = f dz b , (f/f p0 yz) ) (3) 
yields b f (2f p ) = 4 x 10“ 10 rad 2 /m at 1 AU, This agrees 
very well with the values 6'(2/ p ) ~ 2 — 6 x 10“ 10 rad 2 /m 
observed at 1 AU [e.g., Lacombe et ai, 1988]. 

Cautious extrapolation of this model to the outer 
heliosphere, for the purpose of making first estimates 
of $ sc alone (since P n ( k) is undoubtedly not Gaus- 
sian), is justified by the model’s simplicity and the good 
agreement with observations at 1AU. At 50 AU, with 
h = r p ~ 10 3 km and An/n = 10 -2 again, the model 
predicts b f (2f p ) = 2 x 10 _u rad 2 /m. Then, assum- 
ing b f (f/f p ) is constant along a 50 AU path length 
(the minimum likely distance noted above), (3) predicts 


0 8C (f = 2kHz) = 2.4 rad and 0 sc {f = 3kHz) = 1.0 rad, 
corresponding to radiation near 4 and 6 times the local 
fpOy respectively. For reference, a source with plausible 
dimensions ~ 1 AU in the outer heliosphere [e.g., Cairns 
et ai, 1992] subtends an intrinsic half-angle 0j n ~ 0.01 
rad at a distance of 50 AU. Table 1 compares the scat- 
tering results with Gurnett et ai \ s [1993] data. Two 
results follow. First, the apparent source sizes observed 
are much larger than plausible 0 in values, suggesting 
that an effect such as scattering dominates the data. 
Second, the surprisingly good agreement between model 
and data makes it plausible that scattering contributes 
majorly to the observed source size. 

Angular broadening measurements of astrophysical 
sources along lines of sight within 0.3 AU of the Sun 
yield a second model [Armstrong and Woo, 1980; Span - 
gler and Armstrong, 1990]: the 3-D wavenumber density 
spectrum is Kolmogorov with a normalization Cjy that 
depends on heliocentric distance P; i.e., P„(k,P) = 
Clf(R)k~ a with a = 11/3 and 

/ p \ —4.05 

C%{R) = 3.9xl0 15 (— J m- 20/3 (4) 

Here Rs is the Sun’s radius. Since the solar wind speed 
v sw is independent of R and the background plasma 
density falls off as P“ 2 in steady-state, this model corre- 
sponds to the relative density fluctuations being essen- 
tially independent of R. For power-law density spectra 
the angular distribution of the scattered radiation is no 
longer Gaussian [e.g., Lee and Jokipii , 1975], yet 0 SC can 
be calculated as the equivalent 1/e angular half-width of 
the radiation. Generalizing previous PWE derivations 
to include finite f p q / /, the result is 

02 / \ _ -2 ( 422 f j Cn ( z ) ^ a 

e, c (z) - 7T (4^ J dz (1 _ /po(2)2//2) 

-a/2) 2 2 -“ \ 2/ ("- 2) 
a/2) (a -2)/ 

where T(p) denotes the incomplete Gamma function 
and r e is the classical electron radius. Inserting (4) 
into (5), the integrand is proportional to P -4 05 for 
/ ^ fp\ fch e rate of scattering decreases with in- 
creasing R for a given /. At 1 AU (5) predicts scattered 
source sizes in good accord with the ISEE observations 
quoted above. Table 1 shows the results of extrapolat- 
ing this model into the outer heliosphere using the same 
path length and observation position as for the Gaus- 
sian model. Once again the calculated effects of scat- 
tering are very substantial. In fact, 0 SC is predicted to 
be much larger than Gurnett et ai ’s [1993] observations 
indicate. This argues that scattering is plausibly very 
important, but also suggests that the level of density 
turbulence is smaller than in (4). The differences be- 
tween the Gaussian and Kolmogorov predictions show 
the importance of properly modelling the turbulence. 
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Table 1 . Predicted and observed 0, C (5OAU) in radians 


Model 

9 sc 

9s c 

Observed 0 ac 


(2 kHz) 

(3 kHz) 

(3 kHz) 

Gaussian 

2.4 

1.0 

1.0 

Kolmogorov 

76 

30 

1.0 


Predicted radial variations 


Motivated by the semi-quantitative success of the 
above calculations, and the observation of turbulent 
density spectra that are close to Kolmogorov at dis- 
tances ~ 0.3—10 AU, it is appropriate to provide 
more detailed predictions for 9 8C (R , /, f p o) using a Kol- 
mogorov density spectrum with level Cjj(R) = 3.9 x 
lO 15 ^#/ i?s)" 4 that is closely related to (4). The 
quantity w describes the relative value of Cjf(R) rel- 
ative to (4). The path length integral (5) is easily con- 
verted into an integral over R for a radial path from 
the source to the observer. Writing f p (R) = fp\R~ l in 
steady-state, where f p \ is the plasma frequency at 1 AU 
and R is the heliocentric distance measured in AU, the 
path length integral can be evaluated analytically: 


0*/ 3 = B A n/3 w 


' 2 
A^R 


+ 33 l°8( 


R- A 
R + A 



( 6 ) 


Here the observer is at Ro , the source is at R*, A = 
/pi//, and B — 7.0 x 10 " 12 . It may then be shown that 

0 app OC Ro 1 8 r 2 2 w 0 - 6 (7) 


when Rs » Ro, f ~ 2/ p0 (z), and 6 SC » 6 m . Figure 1 
shows the variations in 0 app (R) predicted by (6) at / = 
2 and 3 kHz for w = 1.0 and 0.1. Note the agreement 
with (7) at large R . At small R ~ 2/ p i//, however, 
9 3C increases faster than Rq 1 8 due to the 1 — f p o(z)/ f 2 
term in (5). Comparing Figure 1 with Table 1 it is 




Figure 2. Comparison of 9 sc (Ro) versus Ro with data. 

clear that scattering is predicted to cause major angular 
broadening of the 2-3 kHz radiation source. 

If scattering is important and the source is in the 
outer heliosphere, then one important prediction of (3)- 
(7) is that flee will decrease with increasing R. Contrar- 
ily, for radiation produced in the inner heliosphere at 
e sc (R) will increase with increasing (decreasing) R 
for R > Ri (R < R\). If scattering is important, then, 
radial trends in 0 app (il) can establish if the radiation is 
produced in the outer or inner heliosphere. (It may be 
necessary to also consider radial variations in 0, n , which 
decreases as the source becomes more distant.) 

These predictions are, of course, preliminary due to 
their reliance on a turbulence model extrapolated from 
within 1 AU and the inclusion of scattering in the so- 
lar wind alone. Future refinements should include new 
turbulence models and (most likely) scattering in the 
heliosheath. Observational tests of predictions for 0 app 
should then constrain the characteristics of the density 
turbulence between the source and observer and, ide- 
ally, the location and nature of the source. 


First evidence for radial trends in the 
roll modulation data? 

Two sets of published roll modulation data exist for 
radiation near 3 kHz and apparent source sizes can be 
extracted therefrom. Voyager 1 was at R ~ 18 AU for 
Kurth’s [1988] data and at R ~ 45 — 48 AU for Gur- 
nett et al } s [1993] data. Kurth’s [1988] data show a 
very weak modulation pattern (the peak and trough 
values differ by ~ 1%) whose asymmetry casts some 
doubt on its interpretation. It is thus appropriate to 
take 0 app = t r rad as an approximate lower limit. In 
comparison, near 48 AU Gurnett et al.’s [1993] Figure 
2 implies significantly larger modulation and therefore 
smaller 0 app . Taking the source to have constant inten- 
sity per unit solid angle, standard calculations imply 
that an upper limit on 0 app is ~ 1.0 rad. (The source’s 





elevation above the spacecraft’s roll plane affects 0 app , 
with larger elevations implying smaller 0 app for the same 
modulation pattern. Here the source is assumed to lie 
in the roll plane.) The two data sets thus show 0 app 
decreasing by a factor ~ 3 from 18 AU to 48 AU. As- 
suming small intrinsic source variations, as suggested by 
the close similarities between the 1983/84 and 1992/93 
events [e.g., Gurnett et a/., 1993], the data imply that 
#app(-R) decreases with increasing ft!, as predicted by (6) 
- (7) and Figure 1. Recent data at R ~ 50 AU [Gurnett 
and Kurth, 1994] have even larger modulations than 
the 48 AU data, thereby further suggesting that 0 app 
increases with R . The lack of roll modulation near 2 
kHz is not inconsistent with scattering theory. 

Figure 2 compares the prediction (6) for / = 3 kHz 
and tv = 10" 2 and 10" 1 with the observed 0 app . Rea- 
sonable agreement is apparent for w = 10" 2 . This im- 
plies three results. First, scattering is indeed quantita- 
tively important in determining the observable charac- 
teristics of the 2-3 kHz radiation. Second, the observed 
variations in 0 a pp(^) are consistent with the source be- 
ing in the outer heliosphere but inconsistent with the 
source being inner heliospheric (e.g., Jupiter). (See 
Gurnett and Kurth [1994] for other arguments against 
a Jovian source.) Third, either the density turbulence 
is weaker by a factor ~ 10 2 than the model (4) or else 
other aspects of the turbulence model or the distance 
to the source need revision in order to obtain semi- 
quantitative agreement between data and theory. In 
either case, these results argue for further work on scat- 
tering theory for the 2-3 kHz radiation and on MHD 
turbulence in the outer heliosphere. 

Conclusions 

Calculations are presented of angular scattering of ra- 
dio waves by density irregularities, with application to 
the 2-3 kHz emissions observed by the Voyager space- 
craft outside ~ 10 AU. Two models for the turbulence, 
both of which are physically plausible and accurately 
account for scattering near and well inside 1 AU, are ex- 
trapolated into the outer heliosphere. For nominal path 
lengths ~ 50 AU and both turbulence models the calcu- 
lations predict significant angular broadening. Indeed, 
scattering alone can account for the source sizes ob- 
served. Preliminary predictions for radial variations in 
the scattered source size 0 9C are presented, based on ex- 
trapolation of the Armstrong- Woo Kolmogorov model 
well beyond 1 AU and a path solely in the solar wind. 
Good agreement between the predictions and data are 
obtained if the turbulence level is ~ 100 times smaller 
than the Armstrong/ Woo model. One prediction is that 
the apparent source size of an outer heliospheric source 
should decrease as R increases. Voyager roll observa- 
tions at R = 18 and 48 AU are qualitatively consistent 
with this prediction, arguing that scattering is impor- 
tant and that the radiation is outer heliospheric and not 
from Jupiter. Additional qualitative evidence for scat- 
tering is provided by the lack of modulation near 2 kHz 
and the variability in apparent source direction. While 
intrinsic changes in the source’s size, direction, and in- 


tensity are possible and even likely, these calculations 
and the analogy with type III bursts etc. demonstrate 
the great plausibility of scattering changing the appar- 
ent size and other characteristics of the source. Scat- 
tering should be carefully considered when interpreting 
and modeling the 2-3 kHz radiation. 
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